Dipolar coupling and multidomain states in perpendicularly polarized nanostructures 
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Stripe states in multilayer systems with perpendicular polarization are investigated by analytical 
calculations within a general continuum approach, applicable to ferromagnetic, ferroelectric, or 
ferroelastic nanoscale superlattices. The competition between the long-range depolarization effect 
and short-range interlayer couplings can stabilize monodomain states and unusual stripe phase 
ground states with antiparallel polarization in adjacent layers. Geometric parameters of stable 
stripe domain states and the phase transitions lines between single domain states, aligned and 
antiahgned stripe states have been derived. The theory is applied to analyze multidomain states 
and phase transitions in antiferromagnetically coupled multilayers [CoPt]/Ru. 
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Multidomain states considerably influence physical 
properties of condense matter systems with spontaneous 
polarization. Such spatially inhomogeneous patterns 
form ground states of ferromagnetic 111, l2|, ferroelectric 
[Si], or ferroelastic Q films. Recently multidomain struc- 
tures have been observed in nanoscale magnetic films 
and multilayers with strong perpendicular anisotropy 
1^. jisl. [T^ . [TtI and in ferroelectric superlattices 
1713. Illl|. Similar spatially modulated states can also 
arise in polar or magnetic liquid crystals [l^ [13], po- 
lar multiblock copolymer layers [llj, in superconduct- 
ing films or magnetic-superconductor hybrids [55] , and in 
shape memory alloy films [23i]. Multilayer systems with 
perpendicular polarization components provide ideal ex- 
perimental models to investigate fundamental aspects of 
ordered structures and stable pattern formation in con- 
fining geometries. Control of such regular depolarization 
patterns is also of practical interest. E.g., patterns may 
provide templates, which could be decorated by nanopar- 
ticles or macromolecules, or they may be used to calibrate 
imaging techniques. Particular interest in ferroelectric 
superlattices is driven by the exciting possibility of using 
ferroelectric nanostructures in nonvolatile memory appli- 
cations, new microelectromechanical systems (MEMS), 
and for nonlinear optics devices 3]. In nanomagnetism, 
antiferromagnetically coupled sup erlattices with strong 
perpendicular anisotropy [l^, 1 14 [T^ are considered as 
promising candidates for nonvolatile magnetic recording 
media and other applications [3]. According to recent 
experiments [sl. [l3l.Tl^ these nanoscale superlattices are 
characterized by new multidomain states, unusual depo- 
larization processes, and other specific effects which have 
no counterpart in other classes of media with perpendic- 
ular polarization. 



'Permanent address: Donetsk Institute for Physics and 
Technology, 340114 Donetsk, Ukraine ; Electronic address: 
a. bogdanov@ifw-dresden.de 

t Corresponding author ; Electronic address: 

lu.roessler@ifw-dresden.del 



We present here a detailed analysis of multidomain 
states in magnetic nanoscale superlattices with perpen- 
dicular polarization components. We derive simplified 
micromagnetic equations for equilibrium parameters of 
stripe domains in nanoscale multilayer systems. These 
mathematical tools provide a clear description of the 
multidomain processes and reveal the physical mecha- 
nism underlying their unconventional properties. It was 
found that in contrast to other bulk and nanomagnetic 
systems, the magnetic states here are determined by a 
close competition between interlayer exchange and dipo- 
lar couplings. The enhanced stray field couplings are re- 
sponsible for the unusual switching processes and specific 
transformation of the domain patterns observed in syn- 
thetic meta mag netic multilayer systems, as [CoPt]/Ru 
and others [U[l3i [3- The depolarization effects re- 
vealed in this paper have a universal character. They 
arise similarly in ferroelectric superlattices and in other 
nanosized polarized media. 

For definiteness we consider stripe domains in mag- 
netic nanolayers. We analyse superlattices consisting 
of N identical layers with thickness h separated by 
spacer layers of thickness s (Fig. [T]). The perpendic- 
ular anisotropy fixes the easy magnetization direction. 
Within the single layers the magnetization Mi(r), may 
be spatially inhomogeneous. The energy can be written 
in a phenomenological approach as 
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where integrals are over the volume Vi of the single lay- 
ers. H'^'^) and Hdi^) are the externally applied and the 
depolarizing magnetic fields, respectively. The unity vec- 
tor n designates the normal to the film. K > is 
a perpendicular anisotropy. The short-range interlayer 
coupling parameter between (nearly) homogeneous lay- 
ers J = J J j{'ri,rj)dvi dvj > (<) favours (anti)parallel 
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orientation of the neighbouring layers. To investigate 
general effects of competing stray field and exchange in- 
ter layer interactions, we consider a simple model of a 
multidomain structure, namely stripe domains with an- 
tiparallcl magnetization of magnitude M = jM^ — const 
within all the magnetic layers. The adjacent domains 
are separated by thin domain walls with energy density 
tj. The magnetic field H is applied perpendicular to the 
layer surfaces. The stripe structure is described by the 
widths of domains d± that are polarized in the directions 
parallel (+) and antiparallcl {—) to the field. The geom- 
etry of stripes in a multilayer structure is conveniently 
defined by the stripe period D = + d-, and a set of 
reduced parameters 

d+ — d^ ^ h s 
9=^7^ 'P = 27r-, v= T = l + v, (2) 
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where q is proportional to the average magnetization in 
a stripe structure, p is the reduced thickness of single fer- 
romagnetic layers, is the thickness ratio between mag- 
netized layer and interlayer, and r fixes the superlattice 
period (Fig. [T]). 

The reduced energy wn = W^Ar/(27rM^7V) of a system 
with N can be written 
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In Eq. ((HI) , the term linear in p is the energy of the domain 
walls. The next two terms are Zeeman energy and the 
short-range interlayer coupling, respectively. The stray 
field energy Wm {p, q) must be derived by solving the cor- 
responding magnetostatic problem Q. The reduced en- 
ergy ^ depends on the two dimensionless materials pa- 
rameters 
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The strengths of the exchange coupling is measured by 
>f given by the ratio between the exchange and stray 
field energies. The parameter A characterizes the bal- 
ance between the domain wall energy and stray field en- 
ergies. It is fixed by the ratio of characteristic length, 
I = a/{4TT M'^), which is a fundamental material param- 
eter, see Ref. I4) and the thickness of the magnetic layers. 

Due to the mathematical identity of electro- and mag- 
netostatic equations [2] the multilayer with stripes can 
be thought of as a set of planes with "charged" stripes 
(Fig. [T]). For two such planes separated by an interlayer 
with thickness a = luH, the magnetostatic energy is 



1 - (-l)"cos(ng) 



exp(— npw). (5) 



The stray-field energy of the "charges" within the same 
plane is = 9^ + /(O). The dipolar coupling energy 
between the "poles" on different sides of the same layer 



■t^ + + + + 

KI^A — 



■ — + + + + + + 




"^-i- + + + + + + + 

+ _+ + + + + + + + + + + ^ 



+ + + + 



FIG. 1: A fragment of a superlattice: two ferromagnetic 
nanolayers (of thickness h) with stripe domains are separated 
by a nonmagnetic spacer of thickness s. The domains are 
coupled by the exchange (J) and magnetostatic forces. 



is /(I). Hence, the magnetostatic energy of individual 
layers is 0, 1] 



wt''^ =q' + /(O) - /(I) 



(6) 



The stray field energy of the multilayer, Wm originally 
derived by Suna,[§] can be written as a sum of "self" 
energies of the layers and interactions between them 



(7) 
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/(Tfc + l) + /(rfc-I)-2/(Tfc). 

(int) 



The upper (lower) signs in Wm correspond to parallel 
(antiparallel) arrangement of the polarization in adjacent 
layers (Figs. [21 [3|). We denote these modes as ferro (F) 
or antiferro (AF) stripes. The factors Fk equals stray 
field coupling energies between two layers separated by 
distance rfc. They are composed of the four contribu- 
tions of the magnetodipole coupling between pairs of the 
planes bounding the layers. In particular, for two adja- 
cent layers, k = 1 (Fig. [1]), the interactions 1^3 and 
2 <-> 4 yield equal (positive) energy contributions /(t), 
while 1^4 and 2^3 yield negative energy contribu- 
tions, /(r + I) and /(r — 1), correspondingly. 

In common polarized systems with characteristic sizes 
far beyond the nanoscale range the equilibrium domain 
sizes are usually much smaller than the individual layer 
thicknesses, p ^ 1. Numerous observations indicate that, 
as soon as domain sizes approaches the layer thickness, 
coercitivity suppresses the formation of regular multido- 
main patterns This establishes a natural limit for 
domain sizes in classical systems. For <C /i there is 
no effective dipole interaction between different surfaces. 
Hence, in Eq.© for all nonzero uj, one hase f{uj) ^ /(O), 
and the stray field energy of the multilayer is reduced to 
Wm^ [T| . In such decoupled superlattices the multidomain 
states should have similar properties as those in isolated 
layers. 

On the contrary, in perpendicular polarized nanoscale 
films and multilayers the periods of regular multidomain 
patterns is of the same order as their thicknesses or ex- 
ceed these thicknesses [1, [l^, • In such systems dipole 
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interactions between different surfaces have a sizable ef- 
fect. Mathematically, this is seen from slowly converging 
sums of interaction terms between poles far apart and on 
different internal surfaces. For such structures, numerical 
evaluation becomes arduous, and sharpened analytical 
methods are required. To overcome the slow convergence 
in Wjyi, we extend the method introduced in [2^ . With 
the help of the identity J^t^™~^^ exp{—nt)dt = m\/n"^, 
the infinite sums in Eq. JS]) can be transformed into in- 
tegrals on the interval [0, 1] 



fiu)~fi0)^7r'il-q')-2pn{u;) 



where 



n{uj) ^LU^ (1 -t)ln 

"'0 



cos^ {TTq/2) 
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dt . (9) 



sinh" {ujpt/2)_ 
Then the dipolar stray field energy ([7]) can be written 



as 



= l~^n{l)±^J2(^-k)E,irk), (10) 
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Sfe(Tfc) = 2f7(Tfc) - n{Tk + 1) - n{Tk - 1) . (11) 

Minimization of wjv with respect to p and q yields the 
equilibrium parameters of the stripes in the multilayers. 
By the form of the multilayer energy Eq. ^ the widths 
of the stripes arc independent of the interlayer couplings 
K. A complete analysis of these equations in applied field 
will be published elsewhere f28| . Here we investigate the 
ground state of the system in zero applied field, H = 0. 
In this case q = and the parameter p is derived from 
the equation dw^ /dp = 



A: 



1 
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fc)4^)(rfc-f 1)) 
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where s[^^(Tfc + l) 217p(Tfc) - r2p(rfc + 1) - f^p(Tfc - 1), 
and ^p{lu) = dVt / dpiuj) . After elementary transforma- 
tions this function can be written 



tin 



tanh 



dt. 



(13) 



Typical solutions of Eq. (fT^ are presented in Figs. [2] 
and [3] 

Ferro stripes. These solutions exist for any layer thick- 
ness Fig. [21 In the limit of thick layers the period D 
increases, while it tends to infinity as h tends to zero. 
The role of the stray-field couplings changes with the 
ratio of the interlayer thickness to the thickness of the 
polarized layers, v. For the limiting cases of small and 
large thickness of the spacer layer given by the limits 
1/ ^ 1 or ^ 1, respectively, the dependence D{h) 
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FIG. 2: The equilibrium reduced period D/l as a function 
of the reduced layer thickness h/l for ferro (F) stripes in an 
A'' = 10 multilayer for different values of parameter v. Inset 
shows the function D/l {h/l) for an isolated layer 01 and 
the parameters of the minimum point. Dashed line indicates 
D oc vTi fit corresponding to Kittel theory. 
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FIG. 3: The solutions D/l {h/l) for ferro and antiferro (AF) 
stripes in systems with u = 0.1. AF stripes exist for thickness 
larger than the critical thickness her- For a two-layer systems 
{N = 2) the stripe period goes to infinity at the critical thick- 
ness. For systems with N > 2 the critical stripe period Dcr 
has finite values. The thin (red) line shows the dependence 
D{h)/l for the period of F stripes in the case A'^ = 4. The 
corresponding solutions p{A) of Eq. (|12p for these parameters 
are plotted in the Inset: F stripes (red) and AF (blue) lines. 
Dashed line indicates unstable solutions. 



approaches a behaviour of isolated layers with thickness 
hN and h correspondingly. In the limit of small ratios 
the stripe period D{h) approaches the period for a 
F stripe state in a single layer with an effective total 
polarized layer thickness hN . For very large separation 
between the polarized layers, i.e., for large ratios v the 
stripe period is determined by the properties of the de- 
coupled single layers with thickness h. The limiting so- 
lution for an isolated layer D{h) (Inset Fig. [J) obtained 
by (see also Ref. 0) has a minimum point with pa- 
rameters {hrmn/l = 0.96067, D^„Jl ^ 16.3136). The 
nonmonotonic behaviour of D{h) reflects the antagonis- 
tic role of magnetic charges in the formation of the equi- 



FIG. 4: The reduced values of the factor g at as a function 
of iV and different values of v. 1 (A), 5 (□), 0.1 (Q)- Inset 
gives functions gN{v) for different values A'^. 



librium stripes. In the case of small domains D h, 
which is typical for classical systems, the dipole interac- 
tion between different surfaces of the layer is negligibly 
small, and only the interaction between charges on the 
same surface give a contribution to the stray field energy. 
For stripes with sizes D > h the interaction between 
charges from different surfaces becomes a noticeable ef- 
fect and counteracts the interactions between charges on 
the same surface. This can be understood as a screen- 
ing effect. As the layer thickness decreases this screening 
effect becomes stronger and suppresses the stray field en- 
ergy. As a result, for h < hmin the extension of domains 
decreases both the sum of domain wall energies and the 
stray field energy, and the domain period increases ex- 
ponentially with decreasing layer thickness. These sim- 
ple energetical arguments demonstrate that in nanoscale 
multilayers, where domain sizes usually considerably ex- 
ceed the thicknesses of magnetic and interlayer, h and 
s, the interaction between magnetic poles through the 
stack strongly influences the equilibrium magnetic states. 
Due to this effect a non-monotonous dependence of stripe 
periods D{h) arises in multilayers between the limiting 
cases of large/small values both of v and h, as shown in 

m 

AF stripes. In this case the solutions of Eq. p(A) 
exist on the finite interval A < Acr and consist of two 
branches with stable and unstable solutions (Fig. [3l In- 
set). Correspondingly the equilibrium stripe states exist 
only in a finite range of the thickness h < her (Fig- El- 
At a critical thickness her the period reaches a critical 
value Dcr for multilayer systems with N > 2 and tends 
to infinity for two-layer system, N = 2, (Fig. [3]). 

In multilayers with a ferromagnetic coupling (J < 0) 
both the stray fields and exchange interactions favour 
parallel arrangement of the magnetization across the 
stack. Hence, the F-stripe mode is the ground state in 
such multilayers. On the contrary, in systems with anti- 
ferromagnetic interlayer couplings, J > 0, the dipole and 



FIG. 5: The phase diagram of magnetic states in variables h/l 
and the reduced coupling x favouring antiparallel orientation 
of adjacent layers for A'^ = 4. Thick lines indicate the first- 
order transitions between different phases. They meet in a 
triple point. The dashed line shows the stability limit of the 
AF stripes. Inset shows a difference between stripe periods 
in F and AF modes. Between the triple point and the critical 
point AF domains are metastable. 



exchange forces have competing character. As a result, 
three different equilibrium phases can exist in such mul- 
tilayers in zero field, see phase diagram Fig. [H For suffi- 
ciently strong interlayer couplings a monodomain phase 
with antiparallel correlations of adjacent layers across the 
stack, and an AF stripe phase exist. The transition be- 
tween the monodomain and AF stripe states for the two- 
layer system, TV = 2, is continuous, as seen from the 
divergence of D{h) (Fig.[3|). For systems with N > 3 the 
transition is first-order. For weaker interlayer coupling 
m:, the F stripe phase becomes stable. The transition be- 
tween the F stripes and the AF stripes or monodomain 
state is a topological transition, and it is always first- 
order. At a triple point, that weakly depends on the mul- 
tilayer repetitions N, all three phases coexist. Lines of 
first order transitions between the AF monodomain and 
AF stripe phase in dependence on N and layer thickness 
h arc shown in Fig. [S) 

In the limit of large domains {D > I) the expansion of 
the integral jl]) for g = yields fJ(tj) = 3/2 - ln(pw/2) 
and energy Q can be transformed in the form 



WN = 1- 



2p 



N 



>^(l-]^)(14) 



where pjq = pN = 2n{hN)/D, An ^ A/N + gN{v)- 
Then, the energy wn (1141) has the form of the energy for 
an isolated layer (see [23|) with "effective" thickness hN . 
This means that in a system with large domains {D ^ h) 
the multilayer stack behaves as an effectively coupled sin- 
gle layer and the domain period approaches the solution 
for a single layer with a total polarized thickness hN with 
an effective characteristic parameter A. 

The function gN{i^) describes the influence of finite 
spacers on the magnetic properties of the multilayer via 
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the redistribution of the internal "charges" within the 
stack. This function is given by 



gN 



-l^J2{N^k)G{iy)-\nN, 



(15) 



fc=i 



where G{i') — 2g{Tk) — g{Tk + 1) — girk — 1) and g{uj) = 
a;^ln((x'). The dependence of gN{i^) on the number of 
layers N and the ratio v — s/h is shown in Fig. 1151 In 
the limit of small v this function behaves as g^ = aj^iv 
where the factor 



N~l 



GN = 2{l-l/N)lnN -m-^ k\nk 



(16) 
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varies from ln(2) for iV = 2 to unity as N tends to infinity. 
The equilibrium domain period in the limit of large TV is 

(cf. M) 



D = 7r/iA^exp(Ajv - 1/2) 



(17) 



Recently, stripe domains have been investigated in 
Co/Pt multilayers [13, O, [ll]- AF stripes have been 
observed in Co/Pt multilayers [l^]- F stripes have been 
investigated in a set of multilayers [Co(4A) Pt {7A)]x 
(with X from 5 to 160) fis'l. The average domain width 
D/2 in our notation plotted versus the total multilayer 
thickness is close to the line typical for an isolated layer 
(Inset, Fig. [2]), the minimum point (/imm, -Dmm) is in the 
region AT « 20 with a period D about 150 nm. This sys- 
tem has Vf = 7/4—1.75. However, the existing data are 
insufficient for detailed analysis and comparison with a 
theoretical dependencies for D{h) in the Fig. [2l More 
substantial results have been obtained on the investi- 
gation of antiferromagnetically coupled (via Ru) ferro- 
magnetic blocks [[Co(4A) Pt (7A)]x-iCo(4AVRu(9A)]Ar 
(with iV = 2 to 10 and X = 2 to 12), Ref. Q. Strictly 
speaking magnetic properties of such ferromagnetic mul- 
tilayers may strongly differ from those of single lawrs (see 
Fig. [2]). However, according to the results of [15] these 
blocks can be modelled by a single effectively ferromag- 
netic layer with total thickness h — IIX — 7 A. Hence, 
i/(X) = 9/(llX - 7) varies from 0.072 {X =12) to 0.6 
(X^= 2). The stripe periods, D = 260 nm for AT = 8 
[15| . are much larger than the layer thicknesses h rang- 
ing from 1.5 to 12.5 nm. Thus, the approach Eq. 
can be applied to describe these multidomain states. In 
particular, the critical line ht{N) of the first order tran- 
sition between the homogeneous antiferromagnetic state 
and F stripes is derived from the following equation 

- l/N)/2 = exp(-7r//(/iiV) - gNf + 1/2) (18) 

and plotted in the Inset of Fig. [S] It should be stressed 
that despite the fact that energy ^1)1 has been reduced 
to the same functional form as the energy for an isolated 
layer, these two model have different physical properties. 
Namely the function gniv) in Eq. (|14p describes the in- 
fluence of the internal "charges" on the equilibrium states 
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FIG. 6: The phase diagram in variables h/l and N with the 
critical lines between the antiferro single domain and the mul- 
tidomain states. The transitional line between F stripes and 
AF single domain state (Inset, solid points) is in a close agree- 
ment with results obtained for magnetic Co/Pt multilayers in 
[1^ . The line with hollow points indicates the transition in a 
model with zero internal charges (t/jv = 0). 



of stripes. The transition line hto{N) between the anti- 
ferromagnetic homogeneous state and F stripes for an 
isolated layer {gN{v) —0) is plotted for comparison in 
the Inset of Fig. [SI The difference between the two line 
(/ito — ht)/hto = sN/{ttI) depends on the ratio s/l and 
increases with increasing N. 

In conclusion, we have investigated a continuum ap- 
proach for the stripe multidomains that arise as regular 
equilibrium depolarization structure of coupled multilay- 
ers. The zero-field phase diagram displays transitions 
between stripe states where adjacent layers are parallel 
due to depolarization effects and other phases where the 
layers are antiparallel under influence of short-range cou- 
plings. In this AF case, it is possible to stabilize mon- 
odomain states in thin multilayer systems. Furthermore, 
antiferro stripe phases can become stable in systems with 
sufficiently strong interlayer coupling. The transitions 
between these phases are first-order. In fact, within 
this model, the ferro stripe structure is metastable in 
the whole parameter range of the model. Therefore, the 
different multidomain phases can coexist in multilayer 
systems and are transformed into each other by domain 
nucleation and growth processes. Hence, severe hystere- 
sis effects can arise in such multilayer systems subject 
to specific coercivity mechanisms. In particular the AF 
monodomain state can arise in the form of two domains 
that transform into each other by a global reversal of the 
magnetization structure. These domains can coexist out 
of equilibrium as remnant of an antiferro stripe state. 
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